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ABSTRACT 

We investigate the effect of primordial non-Gaussianity on the cross-correlation be- 
tween the CMB anisotropics and the 21cm fluctuations from the epoch of reionization. 
We assume an analytic reionization model and an ionization fraction with /nl induced 
scale dependent bias. We estimate the angular power spectrum of the cross-correlation 
of the CMB and 21 cm. In order to evaluate the detectability, the signal-to-noise (S/N) 
ratio for only a single redshift slice is also calculated for current and future observa- 
tions, such as CMB observations by Planck satellite and 21cm observations by Omnis- 
cope. The obtained S/N ratio is 3.2 (2.8) for /nl = 100 (50) in our fiducial reionization 
model. Our work suggests in the absence of significant foregrounds and systematics, 
the auto-correlations of 21 cm is a better probe of /nl than the cross-correlations 
(as expected since it depends on b 2 ), while the cross-correlations contain only one 
factor of b. Nevertheless, it is interesting to examine the cross-correlations between 
21 cm and CMB, as the signal-to- noise ratio is not negligible and it is more likely we 
can rid ourselves of systematics and foregrounds that are common to both CMB and 
21 cm experiments than completely clean 21 cm of all of the possible foregrounds and 
systematics in large scales. 



1 INTRODUCTION 



The inflationary scenario is strongly supported by th e statistical nat ure of the density fluctuations revealed by recent cosmic 
microwave background (CMB) observations strongly 1 Komatsu 20091 ) . The observe d density fluctu ations have an almost scale 



invariant spectrum and nearly Gaussian statistics as it was predicted by inflation l|Komatsull2011 



More recently, the measurement of the degree of deviation form Gaussianity has attracted significant attention as this 
would help pinpointing the correct model among very many inflationary models. For example, the density fluctuations arise 



rev 

from simple slow-roll infla tionary models with a single scalar field is almost purely Gaussian dGuth fc Pilll982 ; Starobinskv 



1982^_ BOTdeenetaL 19831 ). and the deviation from the Gaussianity would then be unobservably small (Rfik^^alJ 



1993 



Gangui et al.l ll994). On the other hand, several inflationary models such as a single field inflation with non-canonical kinetic 
terms or some multi- field inflation models can genera te primordial non-Gaussianities large enough to be observed by ongoing 
surveys, e.g. Planck ( The Planck Collaboration 20061 ) (For comprehensive review see Bartolo et al. 2004 ). 

Wilkinson M icrowave Aniso tropy Probe (WMAP) puts one of the strongest constraints on the local t ype of the primordial 
non-G aussianity 1 Komatsul2009l ) . which is parameterized by the constant dimensionless parameter Jnl as 1 Komatsu fc Spergel 
200 ll ) 



$(x) = $ G (x) + /nl($ 2 g (cc) - 



(1) 



where $ is Bardeeen's gauge-invariant potential, <£>g is the Gaussian part of the potential and ( } denot es the ensem b le ave rage. 
For example, the p resent constr a ints o n the local type of /nl from WMAP are — 18 < /nl < 80 bv lCurto etafl l|2009l) and 
-36 < /nl < 58 bv lSmidt et ail |2009l ) 

The effect of primordial non-Gaussianity appears not only on CMB fluctuations but also on large scale structure. The 
abundance a nd clusterin g of vi rialize d objects is s ensiti ve to the existence of the primordial non-Gaussianity, as it was first 
discussed by I Dalai et al. 

I |2008l ), and Slosar et al. has shown that competitive constraints can be achieved from large 
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scale structure, —29 < /nl < 70. Hig h redshift galaxy s urvey with large volumes are also expected to be good probes for the 
primordial non-Gaussianity ( Desiacques fc Seliak 201p| ). 

With the ongoing and up coming surveys of 21-cm, such as LOFAR (Earker et alJboioh , MWA (|Lonsdale et al.ll2009r i. 
SKA (|Carilli fcRaw lings 2004) and Omniscope (jTegmark fc Zaldarr iaga 2010), we will soon have a map of the 21-cm emission 
line of a large volume of the Universe. 21 cm line emission comes from the spin-flip transition of neutral hydrogen. 21 cm 
fluctuations depends on the abundance and clustering of ionization sources, it is expected to reveal the EoR by the observation 
of the 21 cm fluctuations. Since the fluctuations of the 21 cm line emission trace the large scale structure, they are also expected 
to be sensitive for primordial non-Gaussianity. 

The virialized objects such as first stars and galaxies are considered as possible ionization sources, therefore the constraint 
on the primordial non-Gaussianity can be obtained through the study on the 21 cm fluctuations i nduce d by ionization sources 
evolved from initial conditions with primordial non-Gaussianity. For example, Joudaki et al.l (2011) have investigated the 
21 cm power spectrum with /nl. They showed that SKA and MWA could measure / n l valu es of order 10 and Omniscope has 
the potential to put much more stringent constraint /nl ~ 1- iTashiro fc Sugivamal (|201ll ) have studied the bubble number 



count with /nl on 21 cm maps. With the imminent release of CMB maps from Planck and looking forward to various 21 cm 
experiments, in this paper, we set out to investigate the effect of primordial non-Gaussianity on the cross-correlation between 
CMB temperature anisotropies and 21 cm fluctuations from the epoch of reionization (EoR) in the analytic reionization 
model. 

The cro ss-correlation of CMB-21 cm is expected to be powerful tools for investigating the evolution of the cosmic 



reionization JAlvarez et ah 20061 : Adshead fc Furlanetto 2008), and is expected to be detected by near future observations 



( Tashiro et al.ll201oF . As the removal of foreground and systematic effects in 21 cm is a challenging task, it is beneficial to 



invest in another method other than the auto-correlations of 21 cm. In particular, employing cross-correlations with another 
independent experiment reduces the possible extra-power in auto-correlations due to systematics or foregroun ds within one 



single experiment. The cross-correlation amplitude strongly depends on the evolution of the ionized fraction (jAlvarez et al 
2006). Since the existence of the non-Gaussianity is also expected to amplify the amplitude, it is important to evaluate the 
effect of primordial non-Gaussianity in order to extract the information about the EoR from the cross-correlation between 
CMB and 21 cm fluctuations. In this paper, we focus on the non-Gaussianity effect on the ioniz ed fraction through calculating 
the scale dependent bias of ionized bubbles due to the primordial non-Gaussianity based on lDalal et al. (2008), as ionized 
bubbles are density peak tracers of the density fluctuations in the EoR. This particular formalism is applicable not only for 
21 cm fluctuations but also for other large scale structures which are also density peak tracers. We also calculate the angular 
power spectrum of 21 cm fluctuations with the primordial non-Gaussianity for comparison. 

The outline of this paper is the following. In Sec. II, we present the analytic representation for the cross-power angular 
spectrum of CMB temperature anisotropies and 21 cm fluctuations. In Sec. Ill, we give the simple analytic model of the EoR 
and the ionized fraction bias. In Sec. IV, we show the results of the cross-power spectrum and discuss the effect of primor- 
dial non-Gaussianity. Section V is devoted to the conclusion. Throughout the paper, we use the concordance cosmological 
parameters for a flat cosmological model, i.e. h — 0.7 (Ho = h X 100 km/s/Mpc), fib = 0.05, fi m = 0.26 and og = 0.8. 



2 21 CM— CMB CROSS-CORRELATION DURING THE EPOCH OF REIONIZATION 
2.1 Cosmological 21 cm signal during the EoR 

The observed differential brightness temperature of the 21 cm line from the redshift z in the direction of n is given as in 



Madau et all (| 19971 ) by 



T B (n, z) = T 2 i(z)x H (z)(l + 5 xn (h, z))(\ + 6(n, z)), 



(2) 



where xu is the hydrogen neutral fraction, S XH is the fluctuation contrast of 1h, and 5 is the baryon density contrast. The 
normalization temperature T2i(z) is written as 



T 2 i = 26x H 1 



TCMB \ 



n b h 2 
0.02 



l + z 
10 



0.3 

Sim 



1/2 



mK, 



(3) 



where Tqmb is the CMB temperature and T s is the spin temperature given by the ratio of the number density of hydrogen 
in the excited state to that of hydrogen i n the ground state. Dur ing the epoch of reionization, the spin temperature becomes 
much larger than the CMB temperature ( Ciardi fc M adau 200^). Therefore we assume (1 — Tqmb /T s ) ~ 1 hereafter. 

Applying the spherical harmonic transformation to Eq. (T2J), we can obtain the multipole moments of the 21 cm fluctuations 
a|^. In the linear order, afj^ is expressed as 

af m (z) = 47r(-i) e T 21 (z)D(z)xn J ^3 [S h + 6 Xak ] F/ m (fc), (4) 

where D(z) is the linear growth factor, and 5k and 5 Xil k are the Fourier components of 5 and S xii , respectively. 
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2.2 CMB Doppler signal from the EoR 

The main contribution of the CMB temperature anisotropy on large scales from the EoR comes from the Doppler effect. The 
CMB Doppler effect in the direction n is given in the linear order by 

rva 

T D (n) = -7cmb / dr/Te~ T n ■ v(n,r]), (5) 
Jo 

where 770 is the conformal time at the present epoch, v is the peculiar velocity of baryons, f is the differential optical depth 
for Thomson scattering T(rj) in conformal time, which is given by f = n^a^a, with the electron number density n e , the scale 
factor a and the cross section of Thomson scattering or- The continuity equation gives the relation between the peculiar 
velocity and the density contrast. 

v k = -i(k/k 2 )8 k , (6) 

where the dot represents the derivative with respect to conformal time. 

Using the spherical harmonics expansion to Eq. ((5} with Eq. ((6}, we obtain the multipole components of the CMB Doppler 
anisotropy in the linear order, 

af m = 4%(-i) e j dr] J ^^s T c^BDTe' T ^ -^-j e (k(rio - 77)). (7) 

2.3 Cross-correlation 

The angular power spectrum of the cross-correlation between the 21 cm fluctuations from the redshift z and CMB anisotropies 
Ce (z) can be obtained f rom the ensemble average of both the multipole components . From Eqs. ((4]) and lf7)l. Ci(z) is expressed 



as I Alvarez et al.ll200a ). 



Ce(z) = -T C mbT 2 i(z)D{z)- / dk dr}' [(1 -Xi(z))P S (k) -Xi(z)P.a]jiK*l° - t/(*))lM*(»J0 - r,')] — (Dfe- T ), (8) 

■ Jo .A. dr) 

where Ps is the matter power spectrum, P x s is the cross-power spectrum between 5 and S x . In order to obtain Eq. ([8]), 

neglecting the effect of helium ionization, we assume Xi = 1 — xh and 8 X = — 8 XH . Eq. (|SJ) tells us that the amplitude of the 

cross-correlation strongly depends on the evolution of the ionization fraction through dr/dr). 



3 THE BIAS OF IONIZED FRACTION FLUCTUATIONS 



In order to calculate Eq. ((8j) , it is required to evaluate P x s which depends on the reionization model. Since we focus on large 
scales in this paper, we assume that S x can be written with S and the constant bias b x as 

5 X = b x S. (9) 

In order to obtain b x , we adopt the Stromgren limit case in lAlvarez et al. (2006), in which the ionization fraction is 
determined by the balance between the emitting rate of the ionizing photon and the recombination rate. For simplicity, we 
do not include the effect of the primordial non-Gaussianity in this section. 

First, we divide space into cells of mass m. The different cells have different density fluctuations S. We suppose that the 
ionized fraction in a cell can be written as a function of S, 

ln[l - Xi (S)] =-C(5)fcoii(m,5), (10) 

where £ is the efficiency factor for the reionization process and / co n (m, S) is the fraction of the collapsed mass to mean mass 
density. According to the extended Press-Schechter theory, / co ii(m, 5) is given by 

6 C -6 



/coii(m, 5) = erf 



(ii) 



where S c is the critical density for collapse, cr(m) is the smoothed dispersion of the density fields with the top-hat window 
function associated with mass m and <r m i n is the smoothed dispersion at mass m m i n which is the minimum mass for collapse. 
From the form of Eq. (|10|l . we obtain Xj ~ 1 in the lim it of C/coii ^> 1, while, in the limit of C/coii < 1, we get Xi — C/coii 



which corresponds to Eq. (1) in lFurlanetto et al.l (|2004f ). The efficiency £ is a free parameter and depends on the reionization 
model. For simplicit y, we assume that the i onized fraction is decided by the balance between the ionizing photon rate and the 
recombination rate ( Barkana &: Loebl 2001 ). The ionizing photon rate JV 7 at the over-density region with 8 is obtained from 

e 7 77*71(1 + 5) 



"/coll ($), 



(12) 



where e 7 is the number of ionizing photon par baryon in a collapsed object, 7/* is the star formation efficiency, and At* is the 
star forming duration. The recombination rate is given by 

N rec = a Cl n 2 (l + 8) 2 , (13) 
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X j 




Figure 1. The redshift evolution of the ionized fraction, we adopt £ = 50 and m m i n corresponding to the virial temperature T v i r = 10 K. 



where a is the recombination coefficient and c; is the clumping factor. Therefore, we can write the ionized fraction as 

ln[l - Xi(S)] = = - / coll . (14) 

In this assumption, the efficiency of the reionization is given by 

Since we are interested in large scales, we consider cells associated with large mass m. The rms density fluctuation in 
such a cell is much smaller than unity. Therefore, we can assume that |<5| <C 1 and a(m) <C 1 in most of cells. Applying the 
Taylor series expansion to Eqs. (|11[) and ([15}, we obtain 

/con «/coii(0)[l + (5-1)5], 



where 



C(5)«C(0)(l-5), 



7T fcollCrmin.D(z) 
and C(0) = e-fTji, /aci At+n. 

With these approximations, the background ionized fraction Xi is given by substituting 5 — to Eq. |10 

In[l-5i] = -C(0)/«»n(0) I 

and we obtain the expression of S x as 

S x 

Comparing Eq. ((9} with Eq. (|20p we can write the bias b x on large scales as 

b x 



^—?±hi(l -Xi)(b - 2)5. 

Xi 



i_^ln(l-Xi)(6-2). 

Xi 



(16) 
(17) 

(18) 

(19) 
(20) 

(21) 



Our i-f 'ionization iii<mI<'1 lias two paramoti >]-s. ^ and m m in- In this paper, we adopt £ = 50 and m m i n corresponding to the 

2/3 /^2o M/2 



and m 
b 2001), 



virial temperature T VYr = 10 4 K ( Barkana &: Loe 

-.7 / T rmv i 



Tlmin = 3.3 X 10 



Z+l 



Mr. 



(22) 



io 4 J \ w J V - 147 , 

Figure [1] shows the evolution of the mean ionized fraction. The optical depth in the model is 0.1 which agrees with the WMAP 
results. 



4 EFFECTS OF THE PRIMORDIAL NON-GAUSSIANITY 

Prim ordial non-Gaussianities modify the abundance, merger history and clustering of dark halos (Dalai ct al hood ; Slosaretal 



2008). These effects will evidently induce the early beginning of the reionization process, even though the effect on the 
reionization optical depth is not very strong. For example, the local type non-Gaussianity with /nl = 100 will only enhance 
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Figure 2. The cross-power spectrum for different /nl (top panel) and the ratio of the cross-power spectrum between non-Gaussian 
and Gaussian cases (bottom panel). On both panels, the solid, dotted and dashed lines are for /nl = 100, /nl = 50 and /nl = 10, 
respectively. For comparison, we plot the angular power spectrum for the Gaussian case as the thin solid line on the top panel. The 
dashed-dotted line represents the noise power spectrum in the combination of Planck and Omniscope. 



the optical depth by 1% ( Crociani et al. 20091 ). Therefore we assume that the redshift evolution of the background ionized 
fraction obtained in the previous section is not modified by the existence of the non-Gaussianity. However, as reionization 
sources are collapsed objects in the fluctuation of the density field, the fluctuations of the ionized fraction can be affected 
stro ngly by primordia l non-Gaussianity. 

Dalai et al.l (|200&t ) have studied t he eff ect of /nl on peak heights of the density fluctuations and the bias of density tracers. 



According to Eq. (9) in Dalai et al. ( 20081 ) . we can write the bias of density fluctuations of ionized fraction with primordial 
non-Gaussianity is given by 



&ng = b x + 2(b x — 1)/nl<5b7 



' 2aD(a)k 2 T{kY ^ 
where T(k) is the transfer function (T(k) ~ 1 on large scales), 5b is the critical density for the bias tracer. In this paper, 
the bias tracer is the ionized cell . Therefore we adopt the critical density for the ionized cell, while the critical density is for 
collapse, (5 C , in I Dalai et al. 

Since the ionized fraction is almost unity in an ionized cell, C/coii 2> 1 in such a cell. Therefore we assume that the 
condition for the ionized cell is £/ co ii > 1. This condition gives the critical density for the ionized cell, 

S B =5 C - V2K(()[a^ in - 0-2 (m)] 1/2 , (24) 

where K(() = erf _1 (l - Eq. (HP corresponds to Eq. (4) in lFurlanetto et al.l |2004l ) and Eq. (5) in ljoudaki et all (|201ll ). 

The primordial non-Gaussianity /nl affect f co n(8). However this effect is enough small that we can neglect the modification of 
/nl on /con (5) ( iJoudaki et al.ll201ll reported that the mean critical density Sb is only perturbed by 4% even for /nl = 100) . 
Now we can write the cross-power spectrum between 8 and 5 X with /nl is given by 

Psx(k) = b NG P s (k). (25) 

Using Eqs. JSJ and (|25[) . we calculate the cross-power spectra between CMB and the 21 cm line from z = 10 for different 
/nls. We plot the angular power spectrum of the cross-correlation in the top panel of Figure [2] We also show the ratio of 
the angular spectra between non-Gaussian and Gaussian cases, /?ng = C^ G /C^ ausslan , in the bottom panel. Due to the scale 
dependent bias introduced by /nl, the higher /nl induces higher cross-correlation on large scales, while the effect of /nl on the 
cross-correlation is s mall on smaller scales (£ > 100) and it does not modify the position and the height of the cross-correlation 
peak. As pointed by Alvarez et al. ( 20061 ) . the peak height of the cross-correlation depends on the evolution of the ionized 
fraction. Therefore, these facts suggest that the spectrum of the cross-correlation on large scale have the potential to give 
the constraint on /nl, while one can derive information on the evolution of the cosmic reionization from the peak height and 
position. 

For the purpose of cross-checks and comparisons, we also calculate the angular power spectrum of 21 cm fluctuations 
Cf. According to Eq. @, Cf can be written as 



2 f'°° 

Cf (*) = T 2 i(z) 2 - / dk [(1 - x^fPsik) - 2xi(z)(l-Xi(z))P a t + Xi(*)P**] jf WW - *(*))]> 
n Jo 



(26) 



where P xx is the power spectrum of the ionized fraction and is assumed to be P xx (k) = b^ G P(k). Note that we have neglected 
the redshift distortion here due to the peculiar velocity of baryons. The top panels of Figure [3] shows Cf 1 at z = 10 for different 
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Figure 3. The angular spectrum of 21 cm fluctuations for different /nl (top panel) and the ratio of the angular spectrum spectrum 
between non-Gaussian and Gaussian cases (bottom panel). We adopt the fiducial reionization model (£ = 50, T v i r = 10 4 K). Types of 
lines are same as in Figure [2] The dashed-dotted line represents the noise power spectrum of Omniscope. 




Figure 4. The cross-power spectrum in the rapid reionization model for different /nl (top panel) and the ratio of the angular spectrum 
spectrum between non-Gaussian and Gaussian cases (bottom panel). Types of lines are same as in Figure [2] The dashed-dotted line 
represents the noise power spectrum in the combination of Planck and Omniscope. 



/nl and the bottom panel represents the ratio of the angular spectra between non-Gaussian and Gaussian cases, i?NG- The 
angular spectru m on large scales is also amplified by non-zero /nl, showing similar behavior in the 21 cm power spectrum 



as discussed in (jJoudaki et al.ll201ll ). Compared with the cross-correlation, auto-correlation has a larger amplitude, but the 
degree of the amplification due to /nl is small on scales £ < 10. This is because the term P xx proportional to 6 NG is partially 
canceled by P x g on these scales. On the other hand, the term P xx dominates P x g on large scales £ ~ 2 quickly, since the bias 
A6ng (6jvg — b x ) is proportional to k~ 2 as shown in Eq. I] 23 [I. As a result, the degree of the amplification due to /nl, 7Z, is 
similar to the case of the cross-correlation on large scales £ ~ 2. 

On the other hand, as we can see from Eq. ([8]), the cross-correlation amplitude includes dr/drj, the amplitude of the 
cross-correlation depends on the efficiency of the reionization process. We evaluate the cross-correlation in the following rapid 
reionization model with £ = 400 and T v - n . — 10 5 K. These parameters are motivated by the scenario in which the sources 
of ionization photons are massive objects like QSOs. Although the evolution of the ionized fraction in this model is more 
rapid than in the fiducial model, the optical depth is r = 0.1 which is same as the fiducial one. We plot the angular power 
spectrum of the cross-correlation for the rapid reionization model in the top panel and the ratio of the angular spectra between 
non-Gaussian and Gaussian cases i n of F igure [H 



As discussed in Alvarez et al 



!|20oi l: lAdsIiead fc Furlanettol mplified by the rapidness 

of the reionization process. However, large values of ( and Tvu give lower 8b in Eq. I|23p . This would suggest that bubble 
creations are not such rare events in this rapid reionization model, thus reducing the power of using bubbles as highly biased 
tracer to differentiate between different frmNL- 
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4.1 The Detectability of the Cross-Correlation 

In this section, we calculate the signal-to-noise (S/N) ratio of the cross-correlations in order to study the detectability of the 

cross-spectrum signal with /nl- 

If we assume that foreground correlation between CMB and 21 cm can be removed 1 McQuinn et al. 20061 ; Morales et al 
2006), the S/N ratio for the cross-correlation can then be expressed as 



= ^ J^.^ + 1} \Cf- a ? + W + V)(Cf MB + Nf MB )' (2T) 

where the superscript 21 stands for 21-cm fluctuations, the superscript CMB stands for the CMB anisotropy, and Ci and 
Ni are the signal and the noise power spectrum, respectively. In Eq. (|27p . / s k y is the sky fraction common to both CMB and 
21 cm observations. In this paper, since we are interested in large scales, we consider Planck as CMB observation whose sky 
fraction is almost unity. Therefore / s k y corresponds to the one of the considered 21 cm observation. In the Planck configuration, 
compared with the CMB signal, the experimental noise is very small on our scales of interest, Cf MB 3> A^ MB . Therefore, we 
neglect Np MB in the calculation. 

The noise power spectrum of 21-cm observations is given by 

iVf = -^f%^^) 2 . (28) 
tobsA^ \ 2tv A e g ) 

where Ai/ is the bandwidth, t b s is the total observation time, and ^ max = 2nD/\ is the maximum multipole associated 
with the length of th e baseline D. The sy stem temperature T ays is dominated by sky temperature which is expressed as 
T sys = 2.7(1 + z) 2 ' 3 jBowman et all 120061 ). Aotai is the total effective area which is assumed to be ^4totai — NA B a with 



N being the number of the antenna and A e g bei ng the effective area of one an tenna. For 21 cm observation, we consider 



an optimistic experiment Omniscope (or FFTT) ( Tegmark fc Zaldarriagal 20091). We e valuate the observational noise with 



N — 10 6 , A e g — 1 m 2 , D = 1 km, t Q bs = 4000 hours and / s k y = 2ty ( Mao et al. 2008t ). We plot the estimated noise power 
spectrum as the dashed-dotted lines in Fig. [2] and [4] 

We also show the noise power spectrum of Omniscope in the auto-correlation of 21 cm fluctuations in Fig. [3] Compared 
with the auto-correlation, the noise of the cross-correlation is large in particular, on small scales, because CMB Doppler signal 
from the EoR is proportional to fc -1 and the primordial CMB signal gives large noise on small scales. 

We then calculate the S/N ratio for the cross-correlation for the fiducial reionization case (£ = 50 and T v i r = 10 4 ). The 
S/N ratio has weak dependence on /nl due to the amplification on large scales. We found that the S/N ratio is larger than 3 
for /nl = 100 and is less than 3 for /nl = 50. This fact suggests that any detection of excess power in cross-correlation with 
relatively high S/N ratio implies the existence of high /nl. 

On the other hand, the Omniscope auto power-spectrum S/N ratio is large and reaches S/N ~ 14 for £ m ax = 20. For 
comparison, we eval uate the S/N ratio in SKA, where we adopt N — 1400, A e g = 45 m 2 , D = 0.01 Km and f s k y = 0.0056 



( Joudaki et al.ll201lf ). The S/N ratio is 3.8 for i max = 50 for /nl = 100, although the S/N ratio is 1.2 for £ max = 20. This 
suggests that in the absence of significant foregrounds and systematics, the auto-correlations of 21 cm is a better probe than 
the cross-correlations (as expected since it depends on b 2 ), while the cross-correlations has only 1 factor of b. Nevertheless, 
it is interesting to look at the cross-correlations, since it is more likely we can rid of systematics and foregrounds that are 
common to both CMB and 21 cm experiments. 

The rapid reionization case with £ = 400 and T v i r = 10 5 K bring high S/N ratio. Our estimated S/N ratio for the rapid 
reionization case is ~ 4 for each /nl shown in Figure [4] However it is difficult to distinguish each model at that detection 
level, because the dependence of the amplitude is small as shown in the bottom figure of Figure [4] 



5 CONCLUSION AND DISCUSSION 

In this paper we have studied the potential of the cross-correlation between CMB temperature anisotropies and 21 cm 
fluctuations from EoR to constrain the primordial non-Gaussianities. Assuming the analytic reionization model, we have 
utilized the effect of primordial non-Gaussianitiy on the ionization fraction. We have calculated the cross-correlation to the 
linear order and shown the angular cross-power spectrum, while contrasting against 21 cm auto power-spectrum. 

Due to the scale dependent nature of the effect of primordial non-Gaussianity, the effect is larger at large scales. The 
higher /nl become, the more the angular power spectrum is enhanced, and the enhancement is more significant in lower 
multipoles. Since the amplitude of the cross-correlation depends on the efficiency of the reionization, we also investigated 
the effect of different reionization models on the cross power-spectrum. The overall amplitude of cross-power in the rapid 
reionization model is higher than the overall amplitude of cross-power in the slow reionization model. 

The degree of the amplification due to /nl in the cross-correlation is larger on I > 10 than the corresponding scale in 
auto-correlation of the 21 cm fluctuations. However, the CMB Doppler signal becomes small on small scales and is dominated 
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by the primordial CMB signal. This makes the noise large and the detection of the cross-correlation difficult, when compared 
against the auto-correlation of the 2f cm fluctuations. 

To access the detectability, we have calculated the signal-to-noise (S/N) ratio of both auto- and cross-power of Omniscope. 
In the case of the fiducial (slow) reionization model, the S/N ratio is 3.2 for /nl = 100 and less than 2.8 for /nl = 50. In 
contrast, the S/N ratio in the rapid case is larger than 4 for all /nl > 0. However, as mentioned above, there is no clear 
dependence on the value of /nl > 0. In comparison, the S/N ratio for the auto-correlation is quite large. It reaches S/N ~ 14 
and, even for SKA, the S/N ratio becomes 3.8. This suggests that in the absence of significant foregrounds and systematics, the 
auto-correlations of 21 cm is a better probe than the cross-correlations (as expected since it depends on 6 2 ), while the cross- 
correlations contains only 1 factor of b. Nevertheless, it is interesting to look at the cross-correlations, since it is more likely we 
can rid of systematics and foregrounds that are common to both CMB and 21 cm experiments than completely clean 21 cm 
of all of the possible foregrounds and systematics in large scales. In the calculation of the S/N ratio, we ignore the foreground 
contamination of the cross-correlation between CMB and 21c m fluctuations. In reality, so me of the foregrounds for 21-cm 
observations also have the correlation with CMB observations l Adshead fc Furlanetto 2008). This may affect significantly the 
detection of the signal. Therefore a better model o f the foreground is essential for any 21 cm constraint on the non-Gaussianity. 
The tidal approach suggested bv lPen et al. (2012) can be a potential technique to reduce such foreground contamination from 
21 cm mapping. 

In the paper, we consider only one redshift slice for the 21 cm observation in this paper. We can observe many redshift 
slices by choosing the observation frequency. According to Eq. ((27]), taking many redshift slices would increase S/N ratio 
for each /nl- As a result, multifrequency observation of 21 cm fluctuations can bring the better constraint on /nl than the 
one redshift slice such as considered in this paper. In particular, the S/N ratio for the auto-correlation receives benefit richly 
from multifrequency obs ervation. It can be ex pected that even the S/N ratio for SKA become enough to measure the non- 
Gaussianity as studied in lJoudaki et al.l (|201lh . On the other hand, the sign al of the cross-correl ation in the redshift evolution 



reaches a peak during the epoch when the ionized fraction becomes a half jAlvarez et al.| [20061. In particular, in contrast to 
the case of the rapid reionization, there is possibility to utilize many redshift slices for 21 cm fluctuations in the case of the 
slow reionization, as the cross-correlation signals arises during a long period due to the slow evolution of the ionized fractions. 

Finally, we focus the signals from the Eo R on large scales. However it is well-know n that the distribution of the bubbles af- 
fect CMB anisotropics and 21 cm fluctuations ( Iliev et "aPboOrJ : Furlanetto et al ]|2004l ). Therefore, the effect of non-Gaussianity 
can be expected to arise on small scales. We will leave this to future work. 

In conclusion, our work suggests that in the absence of significant foregrounds and systematics, the auto-correlations of 
21 cm is a better probe of /nl than the cross-correlations (as expected since it depends on b 2 ), while the cross-correlations 
contains only 1 factor of b. Nevertheless, it is interesting to examine the cross-correlations between 21 cm and CMB, as the 
signal-to-noise ratio is not negligible and it is more likely we can rid ourselves of systematics and foregrounds that are common 
to both CMB and 21 cm experiments than completely clean 21 cm of all of the possible foregrounds and systematics in large 
scales. 
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